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MaGICC Halos: Confronting Simulations with 
Observations of the Circum-Galactic Medium at z—0 
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ABSTRACT 

We explore the circumgalactic medium (CGM) of two simulated star- forming galaxies 
with luminosities L ^ 0.1 and 1L* generated using the smooth particle hydrodynamic 
code GASOLINE. These simulations are part of the Making Galaxies In a Cosmo- 
logical Context (MaGICC) program in which the stellar feedback is tuned to match 
the stellar mass-halo mass relationship. For comparison, each galaxy was also simu- 
lated using a "lower feedback" (LF) model which has strength comparable to other 
implementations in the literature. The "MaGICC feedback" (MF) model has a higher 
incidence of massive stars and an « 2 x higher energy input per supernova. Aside from 
the low mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends 
to approximately the virial radius. A significant fraction of this gas has been heated 
in supernova explosions in the disk and subsequently ejected into the CGM where it 
is predicted to give rise to substantial OVI absorption. The simulations do not yet ad- 
dress the question of what happens to the OVI when the galaxies stop forming stars. 
Our models also predict a reservoir of cool HI clouds that show strong Lya absorption 
to several hundred kpc. Comparing these models to recent surveys with the Hubble 
Space Telescope, we find that only the MF models have sufficient OVI and HI gas in 
the CGM to reproduce the observed distributions. In separate analyses, these same 
MF models also show better agreement with other galaxy observables (e.g. rotation 
curves, surface brightness profiles, and HI gas distribution). We infer that the CGM 
is the dominant reservoir of baryons for galaxy halos. 



1 INTRODUCTION 



In the paradigm of cold dark matter (CDM) cosmol- 
ogy, gravitational collapse leads to the formation of viri- 
alized and bound dark matter halos. It is predicted 
that baryons fall into these halos with the dark mat- 
ter, comprising a mass fraction up to the cosmological 
ratio. Recent analyses indicates that stars comprise less 
than 2 5% of the baryons that should have collapsed into 
halos (IConrov fc Wechslerl l2009l: iMandelbaum et al l [20091 : 
iMoster et alJl2010b iGuo et alJl2010b iMore et alJl201lh . This 
presents a "missing baryon" problem reflected in the un- 
known whereabouts of the baryons, which did not cool to 
form stars. Substantial amounts of gas are found in a diffuse 
and highly ionized medium which is referred to as halo gas 
or the circumgalactic medium (CGM). 



The CGM is known to manifest in at least two phases: 
(i) a warm/hot, collisionally ionized gas with temperature 
near the virial temperature (T >~ 10 6 K); and (ii) a cooler 
phase (perhaps predominantly photoionized) with temper- 
ature T < 10 4 K. The latter phase has been detected in 
21cm emission in our Galax y in the population known as 
the high velocity clouds (e.g. Wa kker &; van Woerden 1997) 
lying at distances of tens of kpc ( Thorn et a l. 2006). The 
hotter phase may be observed via bremsstrahlung emission 
with X-ray telescopes, in the intracluster medium of mas- 
sive galaxy clusters (e.g. lAllen et"aD l200lh , but so far ex- 
tended diffuse emission has only bee n detected around one 
disc galaxy, the massive NGC1961 ([Anderson &; Bregmanl 
1201 lh . Smaller disc galaxies have X-ray bright outflows 
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([Strickland et al. 2004), but halo emission remains unde- 
tected. 

The diffuse nature of the CGM makes positive de- 
tections of emission rare. The majority of empirical con- 
straints, therefore, come from absorption-line analysis of 
galactic halos, alo ng coincident sightl i nes to distant quasars 
and galaxies (e.g. Bowen et al. 1996; Lanzetta et al. 1995; 



A number of groups have now used cosmological simu- 
lations to examine the nature and enrichment of the IGM 
at z ~ and offe r comparisons to qua s ar absorption line 
observations (e.g. ICen &TO striker 2006; Dave et al. 2010; 



IChen Tinker1l2008l : iRubin et af]l20ld : ISteidel et aDl201C 
These observations reveal a cool and frequently metal- 
enriched phase traced by HI Lyman series absorption (e.g. 
IChen fe Lanzettair2 003) and low-ionization metal-line tran- 
sitions (e.g IChen fe Tinker! 120081 : iBarton fe Cookd 120091 : 
IRubin et al.ll2012l ). Importantly, galaxies of essentially all lu- 
minosity and spectral type exhibit signifi cant HI absorption 
out to impact parameters R « 300 kpc ([Wakker Savage! 
2009; IProchaska et al.l l201lh . implying a massive (M ~ 
1O 1O M ) extended CGM. 

A particularly useful tracer of the hot phase in quasar 
absorption line analysis is the OVI doublet, which oc- 
curs in highly ionized and enriched regions of the Uni- 
verse. If collisionally ionized, the gas has temperature T ~ 
10 5 — 10 6 K. OVI is observed along the majority of sight- 
lines through our Galactic halo and is be lieved to trace 
coron al material on scales of 10-100kpc (Sem bach et al.l 
2003). This gas is also associ ated with the CGM of local 
L > .1L* galaxies (IStocke et al.l 120061 : IWakker fe Savage! 
2009). IProchaska et al. I ([201 lh have further demonstrated 
that the extended CGM (R ~ 300 kpc) of L 0.1L* galax- 
ies has a high covering fraction to OVI and can account for 
all of th e OVI detected in the p resent-day universe. Most 
recently, iTumlinson et al.l (|2Q1 lh report a nearly 100% in- 
cidence of strong OVI absorption for the halos of L ~ L* 
star- forming galaxies, with an ionized metal mass that likely 
exceeds that of the galaxies' interstellar media, demonstrat- 
ing that the CGM is a major reservoir of highly ionized 
metals at z ~ 0. 

The observations reveal a multi-phase, highly ion- 
ized and met al- enriched CGM around present-day gala x- 
ies (also see iThom fe Chenl l2008bl : ISavage et all lOllal lbl: 
IWakker et all l2012h . The high degree of enrichment de- 
mands that a large mass of gas (and metals) is transported 
from galaxies and/or their progenitors to the CGM. Thus, 
these observations provide robust constraints on the pro- 
cesses of gas accretion and feedback which may be compared 
directly against models of galaxy formation. 

Previous theoretical work on the CGM for individual 
galaxies has largely focused on analytic or semi-analytic 
treatmen ts of idealized gas and dark matter / temperature 
profiles (iMo fe Miralda-Escudd Il996l: iTelfer et al.l 120021 : 
iMaller fe Bullockll2004l ; iBinnev et alJl2009h . Although these 
works provide crucial insight into the nature of the CGM, 
they lack proper cosmological context (e.g. mergers) and 
have generally not included the role of feedback from star- 
formation in the galaxy and its satellites. The latter is now 
considered critical to match a wide-range of galaxy observ- 
ables including the luminosity function, the stellar mass- 
halo mas s relationship, a nd the baryonic Tully-Fisher rela- 
tionship (|Guo et al.ll201lh . Galactic-scale winds likely result- 
ing from stellar feedback are a ne arly g eneric feature of star 
forming galaxies ([Shaplev et al.l 120031 ; IWeiner et al. 1 120091 : 
IRubin et al.l 12010) and may significantly influence proper- 
ties of the CGM. 



lOppenheimer et al.l 120121 : ISmith et al.l l201ll : ICenl [2012) 
These studies have not explicitly examined the physical na- 
ture of the CGM for individual galaxies and likely have in- 
sufficient resolution to perform such analysis. A few authors 
have examined the distribution of gas within individual 
galactic halos, but have not yet t reated the distribution of 
metals or the impac t of feedback ([Keres &; Hernq uist 2009; 
IStewart et all l201lh . At z ~ 3. iKawata fe Rauchl ([2007 ) 
have analyzed numerical simulations of outflows in L* pro- 
genitors, and found that OVI better refle cts the strength 
of gal actic winds than HI. M ost recently, iFumagalli et al.l 
(|2012h and IShen et al.l (|201lh have examined the metal en- 
richment of the CGM with the latter demonstrating that 
metals can be ejected to large distances at early times from 
L* progenitor galaxies while their gravitational potential re- 
mains low. 

We have recently begun the MaGICC project to use 
sufficient stellar feedback to simulate galaxies that match 
the stellar mass-halo mass relationship. Early results of 
the project have shown that t he ejection of low angu- 
lar momentum gas via outflows (|Brook et al.l 1201 lh redis- 
tribut es angular mome ntum via large scale galactic foun- 
tains ([Brook et al.ll2012h and thus play a crucial role in form- 
ing disc galaxies, particularly those without classical bulges. 
In this way, we form galaxies that match scaling relations be- 
tween rotation velocity, size, luminosity, colour, stellar mass, 
halo mass, HI mass, b aryonic mass and me tallicity (Brook et 
al. 2012, submitted) . iMaccio et al.l ([2012) also showed that 
these galaxies have cored dark matter density profiles. Here, 
we test the baryon cycle of these simulations with the con- 
straints provided by observations of the CGM, and in par- 
ticular the column densities of H I and O VI as observed at 
z = 0. 



2 SIMULATIONS 

We resimulate a suite of simulations drawn from the Mc- 
Master Unbiased Galaxy Simulations (MUGS. IStinson et al.l 
l2010h . The simulations are listed in Table 1. The highest 
mass galaxy, HM, is g5664 from MUGS and is about half the 
Milky Way's halo mass. The low mass galaxies, LM_MF and 
LM_LF, are rescaled versions of MUGS initial conditions, 
each a factor of eight lower mass than the MUGS simula- 
tions, allowin g us to explore ma ss dependence. LM_MF was 
also used in Brook et al. 

The simulations were evolved using the smoothed par- 
ticle hydrodynamics (SPH) code gasoline ( Wad slev et al.l 
2004). Supernova energy is impleme nted using the blast- 
wave formalism (| Stinson et al . 2006). Metals are ejected 
from type II supernovae (SNII), type la supernovae (SNIa), 
and the stellar winds of asymptotic giant branch (AGB) 
stars. Ejected mass and metals are distributed to the 
near est neighbour gas p articles using the smoothing ker- 
nel ([Stinson et al.l l2006l ) . Metal diffusion is included and 
metal cooling is calculated based on the diffused metals 
([Shen et alJl2010j ). 

Two different star formation and feedback models are 
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Table 1. Simulation data 



Name 


M to t a 
(Mo) 


(Mo) 


Luminosity 
(L*) 


c* d 






200 


HM_MF 


7 x 10 11 


1.4 x 10 10 


0.84 


0.1 


0.175 


Q. 





HM_LF 


7 x 10 11 


4.8 x 10 10 


0.79 


0.05 








LM_MF 


1.8 x 10 11 


3.7 x 10 9 


0.15 


0.05 


0.1 




-200 


LM_LF 


8.8 x 10 10 


8.7 x 10 9 


0.13 


0.05 









Low Feedback 



MaGICC Feedback 



a Mtot is the virial mass of the halo including dark and baryonic 
matter. 

^ M* is the total stellar mass. 

c V-band luminosity compared to My = — 21. 

^ Star forming efficiency. 

e Radiation Pressure Efficiency. 

/ Initial Mass Function: C ^Chabrierl (|2003h : K= lKroupa et al.l (|l993h 



employed. The original MUGS simulations formed stars 
when gas reached a density of 1.0 cm -3 , used a Kroupa 
et al. (1993) IMF and deposited 0.4x 10 51 ergs per su- 
pernova explosion. We refer to this as the "lower feed- 
back model" (LF) in this study, but note that this feed- 
back strength is comparable to or even stronger than 
most i mplementations that ar e currently run in the liter- 
ature ([Scannapieco et aT]|2012h . In our "MaGICC feedback 
model" (MF), four changes have been made to our imple- 
mentation of star formation and feedback: 



• we use the more common Chabrier (2003) IMF which 
creates more massive stars for a given stellar mass; 

• the star formation density threshold is increased to 9.3 
cm -3 ; 

• the energy input from supernovae is increased to 10 51 
ergs; 

• we include energy from radiation released by the mas- 
sive young stars before they explode as supernovae. 

Radiation pressure from massive stars can have signifi- 
cant effects on the s cales that are resolv ed in our simulations 
(jNath fe Silkll2009l : [Murray et alJl201lh . Massive stars typi- 
cally produce 10 50 ergs of energy per Mp, yet this couples 
only weakly to the ISM (jFrever et all 2006). To mimic this 
inefficiency, we inject a fraction of the energy as thermal en- 
ergy in th e surrounding gas but do not turn off cooling (for 
details see iBrook et alJl2012h . Such thermal energy injection 
is highly inefficient at the spatial and temporal resoluti on of 
cosmological simulations, and is rapidly radiated away (iKatd 
1992) This feedback is even less efficient at low resolution. 
Following a parameter search designed to match the stellar 
mass-halo mass res olution from halo abundance matching 
(|Moster et alj|2010h , we inject 17.5% of radiation pressure 
to the surrounding gas as thermal energy in the lower reso- 
lution (HM) simulation, but only 10% in the higher resolu- 
tion runs(LM). The overall coupling of energy to the ISM is 
minimal, but sufficient to reduce star formation in the region 
immediately surrounding a recently formed star particle. 

The HM_LF and HM_MF simulations use the same 
initial conditions and each simulated galaxy has an abso- 
lute l/-band magnitude My ~ —20.8 implying a luminos- 
ity L « 0.8L*. This is somewhat surprising given that the 
two runs yield very different stellar masses for the galaxy; 





-400 -200 200 400 
x [kpc] 



-200 200 
x [kpc] 



Figure 1. Column density maps of HI (top) and OVI (bottom) 
of the low (left) and high (right) feedback simulations. The maps 
are all aligned so that the discs of the galaxies are edge on to the 
viewer. The results on scales greater than 10 kpc are relatively 
independent of viewing angle. 



this difference is compensated by the fact that each galaxy 
follows a very different star formation history (see Figure 
[7J) such that the different stellar ages result in comparable 
V-band luminosity, but B — V colours of 0.3 7 for HM_MF 
and .55 for HM_LF. HM_MF was analysed in lMaccio et al.l 
(|2012h . LM_LF and LM_MF use different initial conditions 
of similar halo mass. They are selected because each has 
M v « -19. The LM_LF initial conditions run with MaG- 
ICC feedback have My ~ —16.2. To summarize the internal 
properties of the galaxies, the lower feedback galaxies suffer 
from the problems of angular momentum loss that have long 
plagued galaxy formation simulations: dense central stellar 
bulges, centrally peaked rotation curves, dark matter cusps 
and too many stars relative to halo mass compared to obser- 
vations. The MaGICC feedback simulations result in galax- 
ies which match the stellar mass-halo mass relation, have 
slowly rising rotation curves, and dark matter cores. The 
MaGICC feedback simulations provide significantly better 
matches to the internal properties of observed disc galaxies. 



3 RESULTS 

To study the CGM properties of the simulated galax- 
i es, we find the galaxy u sing the Amiga Halo Finder 
([Knollmann &; Knebd [2009) and define its center and sys- 
temic velocity by the position and velocity of the particle 
with the lowest potential. We may then sample the CGM 
at a range of impact parameters p. To generate surface den- 
sity maps of the CGM, we sum each box over 1 Mpc along 
the line of sight. The maximum velocity of the material in 
the box is \5v\ < 200km s _1 . To estimate the specific col- 
umn densities of HI and OVI, we must estimate the ioniza- 
tion state of the gas along each sight line. Under equilibrium 
conditions, which we assume apply, the ionization state is 
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Figure 2. Temperature-density phase diagrams of the galaxy 
halo simulated with the two different amounts of feedback. The 
top plots are mass weighted by the total gas mass, the middle is 
weighted by the HI mass, and the bottom is weighted by OVI. 
These phase diagrams only include halo gas, that is gas more 
than 3 kpc above or below the midplane or at K X y > 40 kpc. The 
distinct locations of OVI and HI show that the hot and cold gas 
is in different phases in our model. 



determined by the incident radiation field (photoionization) 
and temperature of the gas (collisional ionization). A proper 
handling of radiative transfer effects are beyond the scope 
of this present paper. Instead, we examine the simulations 
in regions expected to correspond to optically thin material, 
i.e. where the HI column density is low because the gas is 
expected to be highly ionized. 

We calculated the ionization states for hydrogen and 
oxygen throughout the high r esolution region as s uming; op- 
tically thin conditions and the lHaardt &; Madaul ([2012) UV 
radiation field evaluated at z = 0. With the Cloudy so ft- 
ware package (vlO.O last described in iFerland et al.1 (Il998h . 
we generated a suite of models varying the density, temper- 
ature, and metallicity of the medium and used the output 
to find the OVI and HI fractions for all the gas in the sim- 
ulation. 

In Figure [T] we present column density maps for HI and 
OVI for the LF and MF runs of galaxy HM. Both simula- 
tions predict a CGM extending to at least 100 kpc traced 
by cool HI gas and the more highly ionized OVI gas. Fig- 
ure [2] presents three phase-diagrams of the CGM material 
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Figure 3. Cumulative gas mass profiles from 500 kpc inwards for 
each of the 4 galaxies. The solid line represents total gas mass, the 
dashed line represents OVI, and the dotted line represent HI. The 
total mass of each component is given in the legend of each plot. 
With the MaGICC feedback, the HI has a more centrally con- 
centrated profile than the OVI. In the ~ L* galaxy, HM_MF, the 
OVI distribution follows the total gas distribution fairly closely, 
but in the lower mass galaxy, the OVI is significantly more ex- 
tended than the total gas, indicating that enriched material is 
blown farther away. 



in each galaxy. The top phase diagrams include the total 
halo gas mass in the simulations. The middle only counts 
HI mass and the bottom only counts the mass of the halo 
OVI. One notes two distinct phases in the top panels: (i) a 
cool (T - 10 4 K), dense (n > 10~ 2 cm 3 ) gas which dom- 
inates the HI absorption and (ii) a warm/hot (T > 10 5 K) 
gas, which the bottom panels show creates the OVI, primar- 
ily through collisional ionization. 

The MaGICC Feedback simulation shows a high den- 
sity hot component, which seems unphysical that accounts 
for 9% of the halo OVI at z = 0. We find that 90% of this 
high density, hot halo gas is a relic of the thermal feed- 
back implementation-the gas has had its cooling temporar- 
ily disabled-and that this gas is concentrated in a ~ 10 kpc 
sphere around the galaxy. This dense, hot gas drives the out- 
flows that populate the CGM with metal enriched material, 
but it represents a negligible fraction of the OVI column. 
We discuss in £13.31 how this dense, hot gas generates a com- 
parable amount of soft X-ray emission to that observed. 

In terms of total mass, Figure [3] shows that the CGM 
of the HM_MF galaxy has 4 x 10 10 M (7 x 10 10 M ) of 
gas to p = 150(300) kpc. (We quote 300 kpc since it 
is the maximum extent at which OVI maintains a 100% 
covering fraction in the simulations.) About half of this 
CGM has T < 10 5 K giving a cool gas mass that ma tches 
very well with previous empirical estimates (|Prochaska et al] 
1201 lh . Regarding metals, the CGM has - 1O 6 M of OVI to 
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R — 150 kpc, again in excellent agreement with the mass es- 
timates for L « L* , star- forming galaxies ([Tumlinson et all 
l201lh . The LM_MF model has a lower CGM mass in both 
gas and metals by about a factor of 3. The LM_LF has even 
less (rsj 10 9 M©) cool gas in the CGM and a negligible mass 
of OVI. This model is a very poor match to the observational 
estimates. 

Returning to Figure [T] one notes that the OVI gas has 
both a smoother and more extended distribution than the 
HI. Again, this reflects the fact that the OVI gas traces a 
hotter and more diffuse phase in the halos of this galaxy. 
In contrast, the discrete clouds of HI are clustered more 
closely to the disc and are the remnant of both gas-rich 
mergers an d cooling out of t h e hot halo. iMaller fc Bullock! 
(2004|) and iKaufmann et alj (|2009h describe how gas can 
cool out of the hot halo. A detailed examination of the evo- 
lution of these cold clouds is beyond the scope of this pa- 
per, but a preliminary investigation suggests that the clouds 
both cool out of the hot halo and are entrained in the out- 
flows. The cooling mechanism may be a numerical artifact 
of SPH whereby particles approach each other due to ran- 
dom excursions and thus create a slightly higher density, 
which makes them cool faster. Fewer condensations are seen 
using hydrodynamics based on th e Eulerian (grid) method 
(|Tevs sier 20021; [Xgertz et al.ll201lh or the new moving m esh 
code AREPO (jSpringellbOlOl : IVogelsberger et alJ l201lh . It 
remains for future simulations using more detailed radiative 
transfer to see if such condensations are an artifact of the 
SPH hydrodynamics scheme. 

In Figure [J we present the surface density profiles of HI 
gas and in Figure [5] we present the surface density profiles 
of OVI gas as a function of the impact parameter p to the 
center of the galaxy. Overplotted on the distributions are 
observed HI a nd OVI column densiti es for sub-L* galaxies 
at z ~ from IProchaska et al.l ([2011) and for sta r- forming 
L* galaxies at z ~ 0.2 from lTumlinson et alj ([201 lh . We also 
mark each panel with the V-band luminosity of the simu- 
lated galaxies and the virial radius, indicated by the vertical 
green line. Considering first the lower feedback HM_LF, we 
find fair agreement for the HI gas although the data are 
systematically higher than the majority of simulated sight- 
lines. The results for OVI are more discordant; the CGM 
of HM_LF underpredicts the observed OVI column densi- 
ties by nearly an order of magnitude at all impact param- 
eters. In contrast, the HM_MF predicts nearly 10 x higher 
OVI surface densities and therefore provides a reasonable 
match to the observations. Furthermore, this model yields 
qualitatively better agreement with the HI observations. In 
these respects, this single simulation of an L « 0.8L* galaxy 
has a CGM with characteristics matching current observa- 
tions. At lower masses, the differences between the high and 
lower feedback models become starker. LM_LF has an unde- 
tectable OVI content beyond the optical extent of the galaxy 
even though the data for galaxies of comparable luminosity 
all show detections of OVI. In contrast, LM_MF predicts 
a factor of 10 4 higher OVI surface densities, matching the 
observations. 



3.1 Halo Oxygen Distribution Function 

Figure [6] shows the oxygen distribution function of the halo 
gas in the four simulations. The MaGICC feedback cases 
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Figure 4. Radial profiles of the column density maps of HI for 
the four simulated galaxies. Large blue do t s are observations of 
0.1 < L < L* from the IProchaska et al.l (|201ll ) galaxy sample 
while large green squares are galaxies with L<0.1L*. The solid 
green line represents the virial radius, r v i r , for each of the four 
halos. 



result in a broad distribution of oxygen abundance with a 
peak at half the solar oxygen abundance compara ble to the 
abundance s measured in Milky Way halo gas ([Gibson et al.l 
2000, 200 ll ). The Lower Feedback simulations are at least 2 
orders of magnitude less oxygen enriched. In total, HM_MF 
has 1.2 xlO 8 Mq of oxygen in its halo, roughly twice the 
7.7 x 10 7 Mq of oxygen contained in the ISM in the disk, 
while LM_MF has roughly an equipartition between the 
1.6 x 10 7 Mq in its h alo and 2.1 x 10 7 M© in its disk (cf 
Tumlinson et "all l201lh . Lower feedback generates gaseous 
halos with significantly less oxygen. 



3.2 Evolution 

To gain a sense of how the CGM developed, Figure [3 shows 
the evolution of total mass, along with stellar, OVI , cold 
(T < 10 5 K), and hot gas (T > 10 5 K) masses within r vir 
for the simulated galaxies. The total halo mass increases by 
two orders of magnitude in the 2 Gyr between z = 6 and 3 
in each simulation. Subsequent merger /accretion events are 
reflected in mass "jumps". With lower feedback, most star 
formation happens early, at the same time the total mass is 
rapidly increasing. High feedback delays this star formation. 

In each galaxy, hot gas develops at the same time as 
star formation, indicating that stellar feedback is the initial 
source of hot gas in galaxy halos. As the halo grows, accret- 
ing cold gas is shocked to the virial temperature, so that by 
z = 0, only one- third of T > 10 5 K gas at z = was directly 
involved in a stellar feedback event. The mass history of the 
OVI follows stellar mass in HM_MF indicating that oxygen 
produced in supernovae readily makes its way into the halo 
as a hot gas. In HM_LF, the OVI mass initially rises, but 
then drops during its major merger at z = 1.5 after which 
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Figure 5. Radial profiles of the column density maps of OVI for 
four galaxies spanning a range of masses simulated using MaG- 
ICC feedback simulations. T he large dots are observ ations of 0.1 
L*<L<L galaxies from the IProchaska et al.l (|201lh (blue) and 
iTumlinson et al.l (|201ll ) (green) galaxy samples. The solid red line 
represents the virial radius, r v i r , for each of the four halos. 




-3 -2 
[O/H] 

Figure 6. The oxygen distribution function for the four galaxies 
in our sample. The MaGICC feedback produces a significantly 
metal-enriched gas halo with a peak in both the HM and LM 
galaxy just above half solar. 



the OVI mass remains constant indicating that the lower 
feedback is unable to drive metal-enriched outflows. 

It is surprising that even though the feedback imple- 
mentation is thermal, the MaGICC feedback does not lead 
to more hot gas, but rather more cold gas. As described ear- 
lier, this is because some cold gas gets entrained in outflows 
and there is rapid cooling of halo gas. We speculate that 
the similar amount of hot gas in HM_LF and HM_MF that 
have different feedback implementations is representative of 
some maximum amount of hot gas that can exist in a halo 
before the gas cools; a larger parameter search is required 
to address this question. 




Figure 7. The mass evolution of the total halo mass (blue), stel- 
lar mass (dotted- green), 10 3 xOVI mass (dotted-magenta), total 
gas (dashed light blue), and hot (T> 10 5 K) gas within r vir for 
the 4 simulations. The thin black line shows the cosmic baryon 
fraction of the total halo mass. 



3.3 X-ray Emission 

One consequence of predicting that the OVI comes from a 
massive halo of collisionally excited gas ranging up to 10 7 
K is that it should produce diffuse soft (0.2-2.0 keV) X-ray 
emission. Observations of disc galaxies have detected X-rays 
from outflows less than 10 kpc from t he disc, but never Lx > 
5.2 x 10 39 erg s" 1 outside 10 kpc dStrickland et al.ir2 004; 
lAnderson Bregmanl 1201 ll ) . Using a naive iNavarro et al.l 
(|l995h estimate of Lx based solely on bremsstrahlung cool- 



ing, HM_MF has L x = 5.2 x 10 39 (4.4 x 10 37 ) erg s" 1 out- 
side 10 (20) kpc in the 0.2-2.0 keV energy range. Even with 
the large amount of mass in the halo, the gas remains too 
diffuse to emit significant X-rays. It is intriguing that the 
X-ray luminosity from inside 10 kpc, Lx = 5.8 x 10 40 , is 
similar to the observed diffuse, soft X-ray emission f or an 
L B = 2.4 x 10 10 L galaxy (e.g. [Strickland et al.ll2004h . but 
a detailed study of this emission is beyond the scope of this 
work. 



4 CONCLUSIONS 

We have studied CGM gas in high resolution galaxy for- 
mation simulations. The stellar feedback is constrained us- 
ing the stellar mass-halo mass relationship. This constraint 
leads to a stellar feedback implementation that decreases 
star formation and ejects hot, metal-enriched gas into the 
CGM. Inside the galaxies, the stellar feedback results in 
slowly rising rotation curves, flattened dark matter density 
profiles and HI mass-Luminosity, baryonic Tully-Fisher rela- 
tions, and mass-metalli city r elations that resemble observed 
galaxies (Bro ok et al.1 l2012h . We compared these galaxy 
models to one with more standard feedback im plementation 
(jStinson et al.1 l2010l : IScannapieco et al.1 l2012h . The lower 
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feedback simulations match none of these commonly ob- 
served internal properties of galaxies. An examination of the 
CGM in our model reveals that: 

• The MaGICC feedback implementation produces ex- 
tended, metal-enriched gaseous coronae that extend to 
impact parameters of ^300 kpc, even around low mass 
(L* ~ 0.1L©) hal os. These results match observed absorp- 
tion l ine features (Wakker &; Savage] 120091 : IProchaska et al.1 
120111 : iTumlinson et al.ll201 lh . 

• The hot coronae extend beyond rw, though most of 
their mass remains contained within r V i r of the simulated 
galaxies. The coronae extend furthest outside r V i r in low 
mass systems. 

• A reasonably strong feedback implementation, com- 
pared with those in the literature (our "lower feedback" 
case) , is insufficient to match HI absorption features at z = 
and severely underestimates OVI , performing particularly 
poorly in low mass (~ 0.1 L*) simulations. 

• The amount of hot gas in the LF and MF runs did not 
differ significantly, but rather: 

- (i) The MaGICC feedback ejected ten times more 
oxygen out of the disc than the lower feedback in the high 
mass case and a factor of 10 4 more in the low mass case. 

- (ii) The MaGICC feedback halo contains ten times 
the cold gas as the lower feedback case in both high and 
low mass simulations. The cold gas results from cooling 
out of the hot halo gas and cold gas entrained in outflows. 

• The MaGICC feedback simulations match the observed 
HI and OVI absorption columns at a range of masses. 

• The total mass of the coronae is several times higher 
than the stellar mass in the simulated MaGICC feedback 
galaxies. Thus, the simulations of individual galaxies pre- 
dict that the "missing baryons" are found in the CGM of 
galaxies. 

In our simulations, we only considered galaxies de- 
fined as star f o rming, by IProchaska et al.l (|201lh and 
ITumlinson et al] ([201 lh . One of the key results of those 
works is the absence of OVI detection in quiescent galax- 
ies. Our simulations are unable to address the question of 
what happens to the OVI when galaxies stop forming stars. 

Our examination focused on the radial distribution of 
HI and OVI gas in the CGM of a simulated L ~ L* 
galaxy and compared these results to observed surface den- 
sity profiles. There are, of course, additional observables 
related to the entire population of OVI absorbers identi- 
fied i n quasar absorption line studies (e.g. iDanforth &; Shulj 
120081 : iThom fe Chenl l2008al : iTripp et al.1 l2008h : any truly 
successful model must also reproduce these resuts. One 
key constraint is the observed distribution of line-widths 
a.k.a. Doppler parameters or b- values) for the OVI gas 
Thorn fe Chenl l2008al : ITripp et al.1 l2008h which gauge, in 
a statistical fashion, the contributions of thermal broaden- 
ing, turbulence, and gas dynamics to the motions of the gas. 
On average, the fr-values for OVI absorbers are substantial 
(b > 20km s" 1 ) and appear to correlate with the strength 
of OVI absorption. This has posed a considerable challenge 
to models where OVI is p redominantly photoionized (e.g. 
iQppenheimer fe D ave 2009), as such gas has a small thermal 
width. Another valuable kinematic constraint is the common 
observation of co-aligned HI and OVI absorption relative 



to a com mon redshift ( approximately half of observed OVI 
systems: IThom &; Chenl l2008bl ). Indeed, this appears to be 
at odds with the apparent separation of strong HI absorp- 
tion (at smaller radii) with the extended OVI absorption of 
our modeled CGM. On the other hand, we do observe non- 
negligible HI absorption to large radii (Figure [4} which may 
trace the OVI gas. A proper comparison of our model with 
these observables will require the generation and detailed 
analysis of line-profiles, as well as the examination of the 
CGM for galaxies spanning a wide range of masses. This will 
be the focus of our next paper, which will consider whether 
the extended CGM of z ~ galaxies can reproduce all of 
the ob served statistics for OVI as argued by IProchaska et all 
(|201lh . 

We further note that the HI surface densities at p < 
50 kpc are sufficiently large that the CGM will yield signifi- 
cant absorption from lower ionization states of heavy metals 
(e.g. Mgll, Sill, Silll). A proper estimate of the column den- 
sities for these ions, however, will require a full treatment of 
radiative transfer (i.e. to account for self-shielding by opti- 
cally thick HI gas). 

The fact that the resultant extended, metal enriched 
CGM around our simulated galaxies matches the observa- 
tions over a wide range of mass provides strong support for 
a vigorous baryon cycle in which outflows and subsequent 
cooling of halo gas play a key role in forming disc galaxies. 
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